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Magnetoelastic properties in field-induced magnetic ordered phases are studied theoretically 
based on a Ginzburg-Landau theory. A critical field for the field-induced ordered phase is 
obtained as a function of temperature and pressure, which determine the phase diagram. It 
is found that magnetic field dependence of elastic constant decreases discontinuously at the 
critical field. He, and that it decreases linearly with field in the ordered phase (-ff > He)- 
We found an Ehrenfest relation between the field dependence of the elastic constant and the 
pressure dependence of critical field. Our theory provides the theoretical form for magnetoelastic 
properties in field- and pressure-induced ordered phases. 
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1. Introduction 

Lattices of interacting quantum spin S = \/2 dimers 
provide a fascinating class of systems to study various as- 
pects of quantum phase transitions.^ For weak interdimer 
interaction, the dimer spins lock into spin singlets and 
have only short range correlations. As these interactions 
grow, a transition to a magnetically long range ordered 
state can appear, depending on the coupling topology. 
The dimer structure allows a useful representation of the 
spin degrees of freedom in terms of dimer basis: one spin 
singlet and three spin triplet states. In the disordered 
phase, the unique ground state is adiabatically connected 
with the simple product of dimer singlets. The elemen- 
tary excitations are gapped magnons and correspond to 
dimer spin triplets. These magnons are hard-core bosons 
which acquire a finite dispersion, because they can prop- 
agate through the lattices by the interdimer coupling. 
The quantum phase transition to the magnetically or- 
dered state coincides with the closing of the excitation 
gap. Hence, the ordering transition can be described as 
a Bose-Einstein condensation (BEC) of magnons.^' ^ 

The quantum phase transition can be induced through 
the variation of some external parameters. One is ex- 
ternal magnetic field which leads to a splitting of the 
spin triplet multiplet by Zeeman coupling. The lowest 
triplet state removes the gap at a critical magnetic field 
and gives rise to an ordered phase. ^' ^ Another parameter 
is pressure which affects the interdimer coupling, in gen- 
eral, strengthening them such that the magnon disper- 
sion increases and the gap eventually closes at a critical 
pressure. 

There are various ways to investigate the phases and 
phase transitions in these spin dimer systems like neutron 
scattering, electron spin resonance, nuclear magnetic res- 
onance (NMR) and magnetization measurements. In ad- 
dition to these, also thermodynamic properties are im- 
portant to understand the phase transition of these sys- 
tems. Ultrasound experiments and the thermal expan- 
sion are the most easily accessible properties and allow 
us to probe the phase transitions. Schmidt et al. have 
measured the field dependence of the elastic constant 



of SrCu2(B03)2 and NH4CUCI3. They found a discon- 
tinuity in the field dependence of the elastic constant at 
the transition with a linear field dependence in the or- 
dered phase. ^^"^"^ They discuss that the field dependence 
of the elastic constants in the ordered phases resembles 
very much that of the uniform differential susceptibility, 
dM/dH. The thermal expansion in TlCuCls, measured 
by Johannsen et al., indicates a strong uniaxial pressure 
dependence of the magnon gap.^"* However, there are few 
theories for the magnetoelastic effect in the field- and 
pressure-induced ordered phases. 

In this paper, we investigate magnetoelastic proper- 
ties of the field- and pressure-induced orders, focusing 
on copper-chrohdes TlCuClg, KCuCla and NH4CUCI3. 
We will introduce a phenomenological formulation for 
the magnetoelastic effects which is the basis of our dis- 
cussion of thermodynamic properties and the Ehrenfest 
relations at the phase transitions. The spin-lattice cou- 
pling relevant for specific materials will be also discussed 
from microscopic point of view. 

2. Ginzburg-Landau theory 

This section is devoted to the formulation of a 
Ginzburg-Landau (GL) theory to describe magnetoelas- 
tic effects in the vicinity of the quantum phase tran- 
sitions. Such theories have been developed in different 
context previously. ^^'^^ We extend the Ginzburg-Landau 
expansion of the free energy to the case of field- and 
pressure-induced (quantum) phase transitions. The or- 
der parameter describing the spontaneously broken sym- 
metry is the staggered magnetic moment m which lies 
perpendicular to the applied magnetic field. We require 
that this theory gives a description of the basic phase di- 
agram in temperature, magnetic field and pressure as we 
know it to some extent for TlCuCla . The free energy den- 
sity for the field- and pressure-induced ordering is given 

by 



F = Am^ 



1 2 
^coe 



(1) 
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where the first and the second terms constitute the mag- 
netic free energy. The latter three terms describe the lat- 
tice with e as the strain which may represent the volume 
or a uniaxial deformation from the equilibrium position 
in the absence of an external pressure P. cq is the corre- 
sponding elasticity constant, and the last term describes 
the thermal expansion of the lattice. Tg is a characteristic 
temperature which can be determined by experiments. 
The coefficient B is a positive constant and A = deter- 
mines the phase boundary: 



ao 



1 



H 



+ 7ie + 72e 



(2) 



with ao > constant. Hq is a characteristic magnetic 
field which can be determined by experiments. As we 
will see later, Hq is the critical field for zero temperature 
and pressure. The field dependence in A is quadratic; to 
lowest order due to time reversal symmetry of the disor- 
dered phase. The coupling to the lattice is given by an 
expansion in e. Note that in order to achieve pressure 
induced order, 71 > is required, since under pressure e 
is negative (e < 0). The exponents a and of T in eqs. 
(1) and (2) have to be chosen to satisfies Nernst's the- 
orem. Since the entropy S = —dF/dT and the thermal 
expansion ai = de/dT should vanish at T = 0, we find 
the conditions: a > 1 and (p > 1. 

2.1 Basic magnetic properties 

The free energy has to be minimized with respect to 
m and e for given temperature T, external field H and 
pressure P. First, we assume for simplicity that the pres- 
sure determines e. We determine the staggered moment 
m and the uniform magnetization M by the following 
equations: 



dF dF 2aom2 H 

= 0, M= — 



dm? 
which lead to 



dH Ho Ho' 



(3) 



m 



2B 
1 

2B 



H'o 



|aoT| + ao 1 TF 
no 



P<Pc, H>H^ 



P>Pc, H>0 



(4) 



and 



r al 2H^5H + 3H^dH'^ + SH^ 



Table I. Exponents for the staggered (m) and uniform (M) mag- 
netizations obtained by the GL theory, mo is the induced stag- 
gered moment at -ff = above the critical pressure. 
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P = Pc 


P>Pc 


m (X {H - 




m cc H 


m — mo oc H^ 
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He 


M (xH^ 


M (xH 



Hc = Hc{T,P)=HoJ^. 

ao 



(6) 



M ~ i 



BHo 

ao 



BH 



Hi 



P<Pc, H>H, 



P>Pc, H>0 
(5) 



where 



doT — ao 



Here, Pc and He represent critical field and pressure, re- 
spectively. Depending on the applied pressure (or strain), 
aoT is cither positive (P < Pc) or negative (P > Pc). The 
first case (P < Pc) leads to field-induce order, when H 
exceeds the critical field He. Here we find M (x H — He 
and m cx {H — Hc)^^^. For the second case (P > Pe), 
order appears for zero magnetic field, giving the follow- 
ing field dependences: M oc H and m — mo oc H"^, where 
"^o = (WotI/'^B)^^^ is the staggered moment at H = 0. 
At the boundary between the two cases, P = Pe, the crit- 
ical field is reduced to zero Hc{T, Pc) = 0, since oot — 
in eq. (6). At the critical pressure, we find M oc H^ 
and m (X H. These field dependences in eqs. (4) and 
(5) (exponents listed in Table I) do not depend on the 
GL parameters {Ho, To, Po, (j) and a), and are consistent 
with the results obtained from the bond-operator theory 
for the interacting spin dimer system.^" Thus, the GL 
theory with free energy (1) is equivalent to the bond- 
operator theory in the vicinity of the phase transitions 
and otherwise at least usefid for a qualitative analysis. 
Very recently, Goto et al. performed magnetization mea- 
surements under pressure and confirmed that the uni- 
form magnetization actually behaves as M ~ H^ at 
P = Pc.^^ We note that the recent Monte Carlo simula- 
tions by Nohadani et al. also support these exponents. 

2.2 Phase diagram 

In this subsection, we demonstrate that the phase dia- 
gram is obtained by giving appropriate values for the GL 
parameters Ho, To, Po and (6. As mentioned above, the 
phase boundary as a function of temperature, field and 
pressure is determined by the vanishing second-order co- 
efficient A. We assume that the system is disordered for 
T = 0, H = and P — 0. Considering magnetoelastic 
aspects of the phase diagram, we neglect for the mo- 
ment the term 726^. Hence, the critical pressure is easily 
obtained by minimizing first the lattice part of the free 
eriergy with respect to e, which leads to 

P 

Co' 

From this, the critical pressure for zero temperature and 
field is obtained as 

aoco 



(7) 



■7ie + 72e , 



Pc{0,0)=PciT = 0,H = 0) 



71 



Po. (8) 



SH = H- He, 
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Here, Pq can be determined by experiments. The critical 
field for zero temperature and pressure is given by 



J?e(0,0) = i?e(T = 0,P = 0)=ifo. 



(9) 



Therefore, Hq introduced in eq. (2) is the critical field 
for T = and P = 0. In the finite (T, P, iJ)-diagram, the 
boundary is determined by 

2 n 



= 1. 



(10) 



Hy p fT 

The critical field, temperature and pressure are expressed 
as 




Tn{H,P)=To 



PciT,H) = Po 



1 + 



H 



(11) 



The parameters Hq, Pq, Tq and </> can be chosen as to 
reproduce the experimentally given phase diagram. In 
Fig.l, the phase boundary lines arc shown as obtained 
from the GL theory with the basic parameters found for 
TlCuCls. The three figures provide cuts for two variables 
fixing the third one. 

At atmospheric pressure (P = 0), the phase boundary 
has the form 



H^{T,0) - Ho = Ho 



'1 + 



(12) 



The exponent (j) can therefore be determined from the 
experimental phase boundary. For TlCuCla, Oosawa et 
al. examined the phase diagram on T — if plane at the at- 
mospheric pressure by specific heat measurements, which 
give the values Hq ~ 5.7 T, ~ 2.1, and To ~ 4.6 
K.^° Note that the exponent satisfies safely the condi- 
tion (p> 1 given above. 

The critical pressure at T = and H = Q is deter- 
mined by Goto et al. by measuring the magnetization 
curve: Pq 0.4 kbar,^^ while inelastic neutron scatter- 
ing experiment under pressure by Riicgg et al. reported 
that Pq ~ 1.0 kbar.^ The latter value was obtained by 
examining the pressure where the excitation gap closes. 
There is discrepancy in the critical pressure between the 
two experiments, and therefore, we leave Pq as a param- 
eter and renormalize the pressure P in Fig. 1. 

The temperature dependence of Hc{T, P) for fixed P is 
shown in Fig. 1(a), where the ordered phase lies above the 
phase boundaries. Figure 1(b) displays the P dependence 
of Hc{T,P) for fixed T, where magnetic order appears 
above the lines. Increasing the temperature enlarges the 
disordered region. Figure 1(c) shows the phase bound- 
aries in the P — T plane for fixed magnetic fields. The 
ordered phase lies below the lines. 
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Fig. 1. (a) T — H phase diagram for fixed pressure. Values of 
pressure are on the figure in unit of the eritieal pressure Pq. (b) 
P — H phase diagram for fixed temperature. Values of tempera- 
ture axe on the figure in unit of Kelvin, (c) P — T phase diagram 
for fixed magnetic field. Values of magnetic field are on the figure 
in unit of Tesla. Parameters are chosen as Hq = 5.7 T, cp = 2.1 
and To = 4.6 K. 



The result of </> = 2.1 is in apparent conflict with 

the theories discussing the imiversal properties of BEG 
which leads to ^bec = 3/2.'^ However, the temperature 
range, where the universal exponents are valid, may be 
rather small as was found from extensive numerical sim- 
ulations by Nohadani et al.^^ Very recently, Kawashima 
confirmed this point by studying an effective XY model 
with numerical simulations.^^ Misguich et al. also exam- 
ined this point by the mean-field theory of hard core bo- 
son of excited triplet magnons with realistic dispersion 
relations. ^"^ They reported that the critical field behaves 
as Hc{T,0) — Ho ^ only at quite low temperatures 
and that the </> = 2.1 can be seen in a wide range of 
temperatures. 
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3. Ehrenfest relations 

In the previous section, we examined the basic mag- 
netic property and phase diagram, and showed that the 
results are consistent with the recent experiments and 
numerical calculations. In this section, we study proper- 
ties of thermal expansion and elastic constant within the 
GL theory and derive Ehrenfest relations satisfied at the 
phase transition. We note that the Ehrenfest relations do 
not depend on details of the parameters Hq, Tq, Pq, <p 
and a introduced for the GL free energy. The Ehrenfest 
relations help us study the magnetoclastic effects which 
take place at the field- and pressure-induced magnetic 
phase transitions. 

3.1 Thermal expansion 

So far, we have discussed the influence of pressure or 
lattice deformation on the phase transition and the mag- 
netic order. The feedback of the magnetic order on the 
lattice has been ignored so far. This aspect of magnetoe- 
lasticity shall bo now included, since it c;an be measured 
with high sensitivity and provides helpful probes to char- 
acterize the system further, even giving access to certain 
microscopic features as we will show below. 

We examine this by minimizing the free energy with 
respect to e in addition to the staggered moment m. We 
obtain the following behavior: 



K 
Co 

K 
Co 



disordered 



ao7i - Hi 
2Bco' 



ordered 



(13) 



Here, we neglect the quadratic coupling term with 72 
which docs not affect the result qualitatively. The lattice 
expands with temperature as described by the first term. 
The strain shows an abrupt change in the temperature 
and field dependence at the onset of magnetic order. For 
a fixed temperature, the magnetoelastic effect appears as 
a linear field dependence in e close to the phase transi- 
tion, since iJ^ — 2± 2Hc{H — He). With increasing 
pressure. He decreases as in eq. (11) and becomes zero 
at the critical pressure for a given temperature. In this 
case, the strain shows a quadratic field dependence, since 
ifc = in eq. (13). 

The first derivative of e in both H and T is discontin- 
uous. We derive the thermal expansion from the strain e 

by 

de 

Oil = 



ST' 



which leads to 






r aK / 






cqTo \ 




ai = < 








aK / 






. CqTo \ 





given conditions on the exponents a, cj) > 1. At the Neel 
temperature T/v, we find a discontinuity in ai which 
we can bring into relation with other thermodynamic 
quantities. For this purpose, we express some of the phe- 
nomenological parameters in the GL free energy by mea- 
surable quantities. The coupling between strain and or- 
der parameter is connected with the pressure dependence 
of the Neel temperature and the critical field. Thus, a 
simple calculation results in 



71 = — 



aoco / dH, 



dP 



aoCo / dT 



N 



rpip 

^0 



dP 



72 = 



2Hi 



Qp2 



ao4 



N 



2T^ 



Qp2 



(16) 



H 



The entropy is given by 
dT 

, dA aKe 



dF dF_dnf_ dFd^ _dF 



dT 



0-1 



(17) 



Here, we used the following conditions to minimize the 
free energy: dF/drrP' = and dF/de = 0. The specific 
heat is obtained by Cp = T{dS/dT)H,p- At the phase 
transition, the first term in eq. (17) leads to a disconti- 
nuity in the specific heat with a value of 

20-1 



ACr, 



,2j,2 



2BTo 



Tn 
To 



> 0. 



This value increases proportional to T'^~^. We may now 
express the ratio a^/B in terms of ACp. 

Now fix the magnetic field H > Hq and consider the 
temperature dependence of the thermal expansion. In the 
ordered phase, it is given by 



(18) 



Oil = 



aK f T 



Tn V dP 



H 



T^J 



(19) 

which goes to zero as T'^~^ for T — > 0. It is obtained by 

replacing the GL parameters by the calculated quantities 
through eqs. (10), (16) and (18). There is the following 
discontinuity of ai at T = T^: 



Aai 



ACr, f dT, 



t-N 



dP 



(20) 



H 



(14) 



disordered 



This connection between discontinuities is one of the 
Ehrenfest relations for the phase transition. 

Analogously, we may consider the situation for fixed 
temperature and increase the magnetic field to cross the 
phase boundary. The corresponding expansion coefficient 
is defined as 



ao7i 



2BcqH^ V dT 



dHt 



012 



ordered 



dH 







ao7i 



BcqHq \Ho 



Hr 



H - He 

Hn 



(15) 

We find that ai vanishes at zero temperature for the 



disordered 
ordered 
(21) 



We note that a2 has a discontinuity at the transition 
{H = He) and exhibits a linear field dependence in the 
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ordered phase {H > He). For varying field, the uniform 
magnetic moment M and the susceptibihty x take roles 
formally equivalent to the entropy and specific heat, re- 
spectively. 



dF , dM 
M = — — — and y — , 



(22) 



where x has a discontinuity at the phase transition: 

Ax=j4f|^V>0. (23) 



Therefore, we find the following discontinuity at the tran- 
sition for a2' 



Aao 



ao7i 



BcqHo Hq 



Ax 



dP 



(24) 



which is an Ehrenfest relation equivalent to eq. (20). 
These Ehrenfest relations can give certain insights into 

the microscopic understanding of the system. In particu- 
lar, the qualitative behavior of Tn and He as a function 
of uniform or uniaxial pressure can be guessed by looking 
at the lattice structure and the involved exchange paths 
as we will discuss below. The Ehrenfest relations can also 
be tested experimentally. Recently, Johanscn et al. found 
a behavior consistent with the given Ehrenfest relations 
in TlCuCla for various uniaxial strains and correspond- 
ing expansion coefficients.^'' 

Very recently, Sawai et al. reported that the uniform 
magnetization and strain have quite similar field depen- 
dence in the field-induced ordered phase. In the field- 
induced ordered phase, we find the following relations in 
the vicinity of the critical field: 



M 



Ax{H - H,), 
K /T " 



OH, 

dP 



{H - H,). (25) 



Here, wc used eqs. (23) and (24). Both Al and e have 
linear field dependence with related coefiicients. This re- 
lation (25) can be tested by experiments. 

3.2 Elastic constant 

Another readily observable quantity is the elastic con- 
stant c which is connected with the sound velocity. At 
the phase transition, this quantity is renormalized and 
shows a discontinuity. We start with the case of field- 
induced order at fixed temperature. In the ordered phase 
{H > He), the staggered moment is given by 



1 

'2B 



ao 



1+ ^ 



H 
H'o 



+ 716 + 726 



(26) 



which includes the magnetoelastic coupling. When we 
substitute this back into the free energy (1), we find cor- 
rection terms of order of which constitute a modifica- 
tion of the elastic constant in the ordered phase. In the 



vicinity of the critical field, it is expressed as 

72^0 H^ - Hi 



C = Co 



^co- 



± 

2B 

± 
2B 



B 



ff2 



272 ao HeH-He 



B Hn 



Hn 



(27) 



The second term gives a discontinuity in the elastic con- 
stant at H = He, which originates from the linear mag- 
netoelastic coupling to the strain (71 em^). The 72 term, 
involving the quadratic coupling to the strain, leads to 
a linear field dependence in the order phase close to the 
phase transition. 

We again replace the unknown coefficients in eq. (27) 
by mcasmablc quantities. Using eqs. (16) and (23), wc 
obtain the following relation in the field-dependence of 
the elastic constant for a given temperature: 



c = Co - Axc^ 



dHe 
dP 



d^Hl 

Qp2 



H -He 



He 



(28) 



Wc note that the discontinuity in the field dependence 
of the elastic constant is negative at the transition. This 
kind of discontinuous change and the linear field de- 
pendence in the ordered phase have been measured in 
NH4CUCI3 for the elastic constant cge by Schmidt et 
^^11-13 pYom their data, they suggested that there is 
a certain scaling of the field dependence of c and the 
differential uniform susceptibility, x = dM/dH: 



2aom2 2aoHdm^ . ( ^ ^ H - He 
X = —EJT- + 0-2 ao- = Ax 1 + 3- 



Hi 



dH 



He 



■ (29) 



However, such a scaling seems only to be attainable un- 
der the following condition: 



d^Hl 

Qp2 



-6 



dHe 
dP 



(30) 



Obviously, this relation is only accidentally satisfied and 
no deeper connection between c and x is found apart 
from the Ehrenfest relation which relates the discontinu- 
ities of c and x at the transition. 

We consider now the case when P > Pq and even at 
zero magnetic field, where order is possible at low enough 
temperature. Then the magnetoelastic coupling results in 
a field dependence of the elastic constant in the ordered 
phase. Using eq. (27), we find 



c{H) - c(0) 



72^0 

B 



H 
H~o 



„2 / a2 tt2 

iix{H)-xm[^ 



(31) 



In this case, the field dependence is more closely con- 
nected with the differential susceptibility. The field de- 
pendence of c allows us also to determine the sign of 72. 

3.3 Ehrenfest relations 

Finally, we turn to the Ehrenfest relations for the elas- 
tic constant. The discontinuity of the elastic constant is 
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Fig. 2. Intradimer and interdimer interactions for TlCuCls and 
KCuCls. a, b, and c axes are along (100), (010), and (001) direc- 
tions, respectively. Positions of dimers are expressed by {x,y,z) 
in the unit of lattice spacing. J (solid line) is a intradimer in- 
teraction. Ji (dashed line) is a interdimer interaction along the 
ladder. J2 (dashed line) and J3 (long dashed lines) are interdimer 
interaction between different ladders. 



(010) 




(0 01) 



Fig. 3. Schematic picture of exchange pathway for intradimer in- 
teraction. 



given by 



Ac= - 



± 
2B' 



(32) 



This can be related to the discontinuities of the expansion 
coefficients in the following way: 
2 

= -Aai 



Ac 



AC„ /dT, 



Tn \ dP 



dP 



<0, (33) 



and 



dP 



-Aao 



dP 



< 0. (34) 



Also these relations can be in principle tested experimen- 
tally. Since Ac < 0, the sign of Aai and Aq:2 is deter- 
mined by the sign of {^^)h and (^^)t, respectively. 

Wc discussed the above Ehrenfest relations for 
TlCuCla which shows no plateau in the magnetization 
curve. In case of NH4CUCI3, it exhibits magnetization 
plateaus at 1/4 and 3/4 of the saturation moment. 
The recent theoretical and experimental studies revealed 
that the plateaus are the consequence of successive phase 
transitions driven by weakly coupled distinct dimers. ^'^ 
Since the each phase transition can be described by the 
GL theory, the Ehrenfest relations (33) and (34) hold 
also for NH4CUCI3. 



3.4 Uniaxial pressure dependence of exchange interac- 
tions 

In this subsection, we would like to link our discussion 
to specific microscopic aspects of the isostructual com- 
poimds TlCuCls, KCuClg, and NH4CUCI3. For this pur- 
pose, we discuss the pressure dependence of the exchange 
interactions identifying the strain in the GL free energy 
with some uniaxial deformations. In particular, wc dis- 
cuss the sign of (^^)_ff and {^p-)T which determine the 
sign of Aai and Aa2 at the transition, respectively. 

The compounds TlCuCl.3 and KCUCI3 consist of two- 
leg ladders which are composed of Cu^+ and Cl^ ions. 
The Cu^+ ions provide localized spin-1/2 degrees of free- 
dom which interact with each other through exchange 
paths via the Cl^ ions. Inelastic neutron scattering ex- 
periments revealed that interladder interactions are sub- 
stantial. ^^"'^^ The system can be modeled by an isotropic 
Heisenberg model with the exchange couplings given in 
Fig. 2. The two-leg ladders, which run along the a axis, 
are located at the corner and center of the unit cell. The 
critical field is a function of intra- and interdimer inter- 
actions'"' 



g^Hc = ^/P - J{Ji + J2 + 2 J3). 



(35) 



Since J and J2 are the dominant exchange interactions 
for TICUCI3, we assume that the pressure dependence of 
the critical field mainly results from the change of J and 
J2. The critical field He is reduced by increase of J2 and 
decrease of J. The Neel temperature Tn is enhanced, 
accordingly. 

In case of the interladder interaction J2, a shortened 
distance between two spins translates into an increase of 
the interaction. J2 is an interaction along (201) direction 
and it would, therefore, increase, if pressure is applied 
along the same direction. In case of the intradimer inter- 
action J, the angle of the super exchange path Cu-Cl-Cu 
is close to 90° (see Fig. 3). The interaction strength more 
strongly depends on these angles than the distance be- 
tween the spins. Application of a pressure along (010) 
direction increases the angle of the exchange path away 
from 90°, and the intradimer interaction J is enhanced 
by the pressure along (010) direction. Prom this, we ex- 
pect 



\ dP 
\ dP 



> 0, 



<0, 



N 



dP 
dP 



< 0, for P II (010) 



H 



>0. for P II (201) (36) 



H 



This determines the sign of the discontinuity of the ther- 
mal expansion ai and 02 given in eqs. (20) and (24). For 
P II (010), Aai < and Aa2 > 0, while they are oppo- 
site for P II (201). These features are consistent with the 
experimental results. 

4. Summary 

Recently, it has been reported that the field-induced 
ordering in TlCuCls shows a weak first order phase tran- 
sition by measuring ultrasonic attenuation^^ and CI site 
NMR."'*'' Measurements of other physical quantities do 
not so far support this finding. While the magnetoelastic 
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coupling, if it were strong enough, could induce a first 
order transition within our theory, we neglect this fact 
here and analyze the Ehrenfest relation for second order 
phase transitions assuming weak coupling. 

In this way, we have examined magnetoelastic effects 
for the field- and pressure- induced ordered phases based 
on the GL free energy introduced by eqs. (1) and (2). 
The obtained exponents for the staggered and uniform 
magnetizations listed in Table I are supported by the 
recent experiments^ and Monte Carlo simulations. By 
giving appropriate values for the GL parameters {Hq, 
To, Po and (j)), we demonstrated that the GL theory re- 
produces the phase diagram obtained by experiments. 
It enables us to have the phase boundary as a function 
of temperature, pressure and magnetic field. Moreover, 
our study describes the behavior of elastic constants and 
expansion coefficients in the ordered phase. Ehrenfest re- 
lations (33) and (34) have been found, relating various 
measurable quantities. The obtained Ehrenfest relations 
do not depend on details of the GL parameters, and will 
help to analyze and understand experiments of magne- 
toelastic effects. A rather simple discussion of the micro- 
scopic structure of spin exchange interaction allows us to 
determine the sign of the discontinuities of the thermal 
expansion for various uniaxial stress via the Ehrenfest re- 
lations. A first comparison with experiments shows that 
the experimentally determined coefficients give a consis- 
tent view within our GL theory. 
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